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Many analytic results for the connectivity, coverage, and capacity of wireless networks have been reported 
for the case where the number of nodes, n, tends to infinity (large-scale networks). The majority of these 
results have not been extended for small or moderate values of n; whereas in many practical networks, 
n is not very large. In this paper, we consider finite (small-scale) wireless sensor networks. We first show 
that previous asymptotic results provide poor approximations for such networks. We provide a set of differ- 
ences between small-scale and large-scale analysis and propose a methodology for analysis of finite sensor 
networks. Furthermore, we consider two models for such networks: unreliable sensor grids, and sensor net- 
works with random node deployment. We provide easily computable expressions for bounds on the coverage 
and connectivity of these networks. With validation from simulations, we show that the derived analytic ex- 
pressions give very good estimates of such quantities for finite sensor networks. Our investigation confirms 
the fact that small-scale networks possesses unique characteristics different from the large-scale counter- 
parts, necessitating the development of a new framework for their analysis and design. 
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1. INTRODUCTION 

In the past, many analytic results on the connectivity, coverage, and capacity of wire- 
less ad-hoc and sensor networks have been obtained. In almost all of the results, it 
is assumed that the number of nodes, n, in the network tends to infinity (large-scale 
networks). In other words, these results are asymptotic. Asymptotic results are very 
important for two reasons. First, they give us good estimates for large-scale networks. 
Second, they show some fundamental trade-offs in the network. However, in many 
practical wireless networks the number of nodes may be limited to a few hundreds 
(small-scale/finite networks). As it is shown in this paper, the asymptotic results cease 
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to be valid for these networks. Thus, it is very crucial from the practical point of view to 
analyze finite networks. These analytic results will essentially help us to understand, 
design, and analyze practical wireless networks, and also to design more suitable com- 
munication protocols. 

For example, consider the capacity analysis of wireless networks 
which has been s tudi ed extensively (e.g., i n (Balakrishan e t al. 2004 



Gupta and Kumar 2000} [Grossglauser and Tse 2001t [Gupta and Kumar 2003 
Li et al. 20011 IPerevalov and Blum 20031 ILiu et al. 200311 ). Today we have a gooc 



a good 

understanding of scaling laws for the capacity of wireless networks. However, suppose 
we need to design a wireless sensor network consisting of an arbitrary deployment 
of a hundred sensor nodes. Some fundamental questions are as follows. What is the 
transport capacity? What are the connectivity and coverage probabilities of such 
networks? How do network parameters such as the communication radius of nodes, 
number of nodes, and so on, affect these properties? Unfortunately, the available 
asymptotic results fail to give answers to these questions. 

To address the aforementioned issues in small-scale networks we need to address 
some inherent problems. First, in large-scale networks we use asymptotic estimates 
that make the analysis much simpler. These estimates are not available in small-scale 
analysis. Thus, small-scale analysis is usually more difficult. Second, even if we can 
perform the small-scale analysis, we usually obtain very complicated formulas that 
are not very useful practically. In this paper, we want to circumvent these problems 
and provide bounds for small scale-analysis. In particular, we are looking for easily 
computable but acceptable estimates for fundamental network quantities. The main 
goal of this paper is to initiate the small-scale analysis of wireless sensor networks. 
To the best of our knowledge, this is the first work to analytically and systematically 
study this issue. 

The main idea is the following. The first key point is to aim at simple and very 
good approximations instead of trying to find complicated exact formulas. To do so, 
we first consider the asymptotic analysis. Some of the estimates in the asymptotic 
analysis are still good for small-scale networks, while others are not. We identify those 
which are not valid and replace them with better estimates. However, this must be 
done carefully, in order to obtain simple and easily computable formulas at the end. 
Specifically, in this paper we list a few important differences between small-scale and 
large-scale analysis. 

As a special case of finite sensor networks, we first study unreliable sensor grids in 
which the sensors are deployed in a grid and each sensor is active with probability p. 
This probability is used to account for both sensor failures and sleeping sensors. A fun- 
damental question is that given an area to be protected, how many sensors should be 
deployed so that every point in the region is covered by at least one sensor (more gen- 
erally, we may require that every point in the region is covered by at least k sensors). 
Equivalently, one can ask if n sensor nodes are deployed in an area, what should be the 
sensing radius of nodes to ensure coverage (or fc-coverage)? The same question can be 
repeated for other network properties such as connectivity and diameter. In this paper, 
we study the behavior of the different parameters in finite sensor grids. We prove that 
all graph theoretic properties of these networks such as connectivity, network diameter 
and capacity, follow a piecewise constant behavior and this is even true for the cover- 
age which is not a graph-theoretic property. This result shows a key difference between 
the behavior of sensor grids and randomly deployed sensor networks, and has some im- 
portant implications from the practical point of view: 1. It shows that increasing the 
communication and sensing radii does not necessarily improve coverage, connectivity 
or any other graph-theoretic property, 2. It suggests that we can completely determine 
the behavior of a vast class of network properties by knowing their values for only a 
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finite number of points. We then find simple lower and upper bounds for the k-coverage 
probability of sensor grids and show that these bounds are adequately close to the real 
value, as an estimate of the coverage probability 

Next, we consider finite sensor networks in which nodes are randomly distributed 
in the unit square. We study fc-connectivity and coverage of these networks. We give 
several results pertaining to these properties. We first show that the previous asymp- 
totic results on coverage and ^-connectivity are not accurate for the finite case. We 
then provide a very simple formula for the ^-connectivity probability of finite sensor 
networks and show that the formula is very precise. We also study the coverage prob- 
ability of random networks where we prove simple lower and upper bounds for the 
coverage probability. 

The remainder of the paper is organized as follows: Next section provides an 
overview of the related work. In Section El we study connectivity and coverage of fi- 
nite sensor grids. Section [3] investigates the fundamental properties of random sensor 
networks such as connectivity and coverage. Finally, Section |4] concludes the paper. 



1.1. Related Work 



Related problems hav e been studied in the context of random graph 
theory [Bollobas 20011, continuum percolation and geometric probability 



[Meester and Roy 1996[ fPenrose 200311, and the study of wireless network 
graphs [Gupta and Kumar 1998] ; IGupta and Kumar 2000] |Xueand Kumar 2_004j 
Booth et al. 20031 IBooth et al. t IShakkottai et al. t IDubhashi etaLl ILi et al. I 



Wan and Yi [ IKumar et al. 20041 . In random graph theory, the model G(n,p) is 



extensively studied, in which edges appear in a graph of n vertices with probability p 
independent of each other. In continuum percolation theory, usually infinite graphs on 
R d are studied. Finally, in geometric probability and the study of graphs of wireless 

netw orks, large-scale graphs over the plane are usually studied. 

In MFranceschetti and Meester 2008[|Freanceschetti and Meester 20 061. the authors 
studied connectivity and critical node life-time for a model of random networks in 
which the density of nodes is kept constant while the area of interest tends to in- 
finity. Furthermore , the throughput scaling of wireless relay networks is studied in 
IDousse et al. 20061 for this model. However, the results in these papers are all based 
on asymptotic analyses and their method can not be applied to the case of finite net- 
works, i.e. networks with finite number of nodes (e.g. less than 1000) on a finite plane. 
In t he analysis of these networks, boundary effects and constant factors (see section 
cannot be neglected as can be for the case of asymptotic analysis. 
The connectivity and fc-connectivity of large-scale wireless networks have 



been investigated in [ Gu pta and Kumar 19981 , BLi et al. 1. IWan and Yi I. 
RPishro-Nik et al. 20041 and IDousse et al. 20021 . In (Dousse and Thi ran 20041 . 



the trade-off between connectivity and capacity of dense networks was ex- 
amined. The transport, information theoretic, and MAC layer capacities 



have extensively been in vestigated (see for examp le IIBalakrish an et al. 2004 



Gupta and Kumar 2000 |Grossglauser and Tse 2001 |Gupta and Kumar 2003 



Li et al. 200H IPerevalov and Blum 20031 ILiu et al. 20031 The grid model for sensor 
networks has also been investigated. In particular, conn ectiv ity, coverage, and d iam- 
eter of sensor grids were studied in IShakkottai et ap . In IKumar et al. 20081 and 
IJanso n 1986], the fc-coverage pro blem for sensor grids and other deployment met hods 
was considered. The authors in [Bal ister et al. 20091 IBalister and Kumar 2 0091 also 
studied coverage for sensor networks in presence of failures and placement errors. 
However, almost all previous analytical results are asymptotic since they consider 
large-scale networks. 
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Analysis of wireless netwo rks with modest number of nodes has generated a lot 
of interest in the recent pas t [Bai et al. 2006; Desai and Manjunath 2002; Gore 2006[ 
IKarmachandani et al. 20061 fYenand Yu 2004"llGhasemi and Nader-Esfahani 20061 . In 
I Desai and Manjunath 2002[ , the authors investigated the problem of connectivity for 
one-dimensional networks (line networks). Using probabilistic method s, they obta ined 
the exact formulation for the probabi lity of connectivity. The auth or of (Gore 2006] pre- 
sented corrections and extensions to [Desai and Manjunath 2002 1. It is noted that both 
of the above cases considered a line network, and the extension to two-dimensional 
networks was achieved by obtaining a loose bou nd using the results from the for- 
mer case. In [G hasemi and Nader-Esfahani 20061 , the authors also consider the line 
network and obtain connectivity results for one-dimensional net works. The threshold 
phenomena for finite wireless networks on a line is studied in BEslami et al. 20101 . 
The authors also find lower and upper bounds on the MAC-layer capacity for such net- 
works. It should be n oted that the main challeng es in finite analysis arise in the two 
dimensional case. In BKarmachandani et al. 20061 . mobility and more realistic models 
were examined. The authors obtained results on t he connectivity fo r both finite and 
asymptotic cases in one-dimensional networks. In HYen and Yu 20041 , some simple lo- 
cal network characteristics such as the link probability (occurrence of a link) and aver- 
age node degrees are studied . The paper also obtains formulas for the average covered 
area. In i Balister et al. 20071 , connectivity and coverage are studied for networks on a 
thin strip of finite length. The authors provide reliable density estimates for achieving 
coverage and connectivity, assuming a Poisson distribution for the nodes. 



2. FUNDAMENTAL PROPERTIES OF FINITE UNRELIABLE SENSOR GRIDS 

In this section we present properties of finite unreliable sensor grids. In particular, 
we prove that a large class of network properties such as connectivity, coverage, and 
capacity can be represented as a piecewise constant function of the communication 
and sensing radii, r t and r s , respectively. We also discuss the implications of this result 
and show the importance of boundary effects in finite networks. We then find an upper 
bound for coverage which can be used to approximate the exact value of the covera ge. 

Here, we consider the sensor network model introduced in IShakkottai et alTl . In 
particular, it is assumed that n sensor nodes are arranged in a grid over a square 
region of unit area. This region is called the deployment region and it is assumed to 
be the unit square centered at the origin. Such a grid is depicted in Fig. [TJ We show 
the deployment region by S - The separation between adjacent nodes is assumed to be 
-7= units. Each sensor node can detect events within some distance from it, called the 

sensing radius r s . Each sensor is active with probability p independently from other 
nodes. The transmission radius of each node is assumed to be r t . In other words, if the 
distance between two sensor nodes u and v is less than r t , then they can communicate 
with each other, thus the edge {u, v} belongs to edges of the graph. It is worth noting 
that our results apply to any deterministic placement of finite sensor networks and 
also any finite deployment region with smooth boundaries. However, for simplicity, 
we consider the above grid model in this paper. We are interested in connectivity and 
coverage. In particular, we assume Pdisc(n,p, r t ) is the probability that the sensor grid 
with parameters n, p, and r t constructs a disconnected graph. We also assume that 
Pcov{n;P,r s ,k) is the probability that each point of the unit square (the deployment 
region) is covered by at least k sensors in the sensor grid with parameters n, p, and 
r s . Thus Pcov(ri,p, r s , 1) is the probability that the whole unit quare is covered by the 
sensor nodes. 
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Fig. 1 . A sensor grid is shown where the active sensors are shown by black nodes, communication and 
sensing radii, rt and r s , are also shown. 

2.1. Sensor Grids: Asymptotic versus Finite Analysis 

We now present some evidence to show that previous asymptotic results diverge sig- 
nificantly from actual values for finite grids. To show this, we consider connectiv- 
ity and coverage. Let us first con sider coverage. The asymptotic coverage probability 
p C ov{n,p, r s ,k) has been found in Kumar et al. 20081. In particular the following fun- 
damental result has been obtained in UKumar et al. 20081 . 

THEOREM 1. (Kumar, Lai, and Balogh 2008) Let e be an arbitrary constant positive 
real number and k be a constant positive integer. Then for n chosen large enough we 
have the following two cases. 

— If r s {n) > \J ^ 1+t ^nf T^~' t nen ^ ne unit square is almost always k-covered completely, 

i.e., p C Qv(n,p ,r s ,k) = 1 - o(l). 
-Ifr s (n) < v / IE5gM 3 then p cov (n,p,r S) k)= o(l). 

Using simulations, authors of UKumar et al. 200 8 1 have shown that this theorem re- 
sults in accurate estimation o{p cov (n,p, r s ,k), when n is large (say n > 10000). Thus the 
theorem is very useful in the design of large-scale sensor networks. Let us now consider 
a sensor grid consisting of 100 unreliable sensor nodes with p = 0.2. If we want to use 
the asymptotic result for this network, choosing e = 0.1, we conclude that if r > 0.229 
then p cov (n,p, r s , k) w 1 and if r < 0.207 then p cov (n,p, r s , k) w 0. We have used ex- 
haustive simulations to obtain an accurate estimate of p CO v{n,p, r s ,k). In Figure |H we 
compare the results obtained by exhaustive simulations and TheoremQ] It is observed 
that the two results differ considerably. For example, at r s = 0.25, the asymptotic result 
predicts that the unit square is covered with probability close to one. However, simula- 
tions show that this probability is only p cov (n = 100, p = 0.2, r = 0.25, k = 1) = 0.018. It 
is clear that for this network the asymptotic analysis cannot provide results that are 
sufficiently accurate. Figure [3] shows that the same situation exists when we consider 
fc-coverage for k > 1. Thus, it is very important to provide finite-size analysis. We also 
observe that the coverage probability obtained b y sim ulations shows several disconti- 
nuities. We prove this phenomenon in the section [2T2l We performed many simulations 
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Sensing Radius r (n=100, p=.2) 



Fig. 2. Comparison of asymptotic results and the simulation results for the coverage probability of a sensor 
grid with p = 0.20. For simulations we picked n = 100 while the asymptotic result hold when n tends to 
infinity. 
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0.4 0.6 0.8 1 

Sensing Radius r (n=100, p=0.2, k=2) 



Fig. 3. Comparison of asymptotic results and the simulation results for 2-coverage probability (k-2) of a 
sensor grid with p = 0.20. For simulations we picked n = 100 while the asymptotic result hold when n tends 
to infinity. 



for different values of n, p, and k to further validate the insufficiency of asymptotic 
results. However, we omit them for brevity 
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2.2. Discontinuity in Properties of Sensor Grids 

Here we prove that a vast class of network properties can be represented by piecewise 
constant functions of r t and r s . We stress that the piecewise property is one of the key 
differences between sensor grids (deterministic deployment) and randomly deployed 
sensor networks. Consider a right-continuous function f(x) : [0,oo] — > BL The function 
/ is said to be piecewise constant if there exists a set of real numbers = x± < x 2 < 
x 3 < and ci, c 2 , such that f(x) = a for all x £ [xi,x i+ i). In this paper we only 
deal with functions for which the number of Xi's is finite. 

Let Q be a property for sensor grids such as coverage, i.e., we say that a grid has 
the property Q if it covers the deployment region. Coverage is an example of geometric 
properties. Another category of properties are graph theoretic properties such as con- 
nectivity. In particular, any sensor grid with parameters n, p, and r t corresponds to a 
graph that can be shown by g(V, E) where V and E are the set of vertices and edges in 
the graph, respectively. The sensor nodes construct V, the set of vertices of the graph. 
There exists an edge between two vertices if their corresponding sensors are within 
the communication range of each other. Any property of g(V, E) is a graph theoretic 
property for the sensor grid. Thus, two different sensor grids will have the same graph 
theoretic properties if they have isomorphic (identical) graphs. We note that coverage 
is not a graph theoretic property. 

Let X = {Xi,X 2 , X k } be a set of points on the plane. Define g(X, r t ) as the graph 
obtained by the following method. The vertices of g are the points in X and there is an 
edge between two vertices X l and Xj, if their distance is less than or equal to r t . We 
prove the following theorem. 

THEOREM 2. Let Q be a graph theoretic property of sensor grids with parameters n, 
p, r t - Let n and p be fixed numbers andpQ(r t ) be the probability that the sensor grid with 
communication radius r t has the property Q. Then pq is a piecewise constant function. 
In particular, there exist = r < n < r 2 < ... < r m < A= f and c\,c 2 , —,c m such that 
PQ( r t) = a ifn g [n--i,n). 

Proof. Let X = {X\,X2, ...,X^} be the set of points in the sensor grid. Let also 
X a c X be the set of active sensors. Assume that g(X, r t ) is the corresponding graph. 
Let p{X a ) be the probability that X a is the set of active sensors, then we have 



It suffices to find = r < ri < r 2 < ... < r m < such that the network graphs 

g(X a , r t ) remain constant for r t £ [r,_i, r,) for any choice of X a and any i £ {1, 2, m}. 
Let D = {dx, da, — , di} be the set of distances between the points in X, and assume that 



< d\ < d 2 < ... < di < -4=. In our grid model we have d. L = * /-, i = 1, I. Then, the 



network graph remains the same when r t £ [di,d i+ i) for any i £ {1,2, ...,m}. This is 
because changing r t within [di, d i+1 ) will not add. or remove any edges. This means that 
we can choose n = di. Thus pq in I© remains constant for r t £ [n-i,n). It is also easy 
to see that pq is right-continuous because the edges in the graphs are formed when the 
distance between two nodes is less than or equal to r t . This completes the proof. □ 

Note that the above discussion shows that any graph theoretic quantity is a piece- 
wise constant function of r t . This includes diameter of the network, MAC layer capacity 




(1) 




(2) 



g(X a ,r t )hasQ 
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BBalakrishan et al. 200411 . fc-connectivity, etc. We now prove that coverage probabilities 
are piecewise constant functions of sensing radius. Note that this cannot be concluded 
from Theorem|2l since coverage is not a graph theoretic property. 

THEOREM 3. Consider a sensor grid with parameters n, p, r s . Let n and p be fixed 
numbers. Then p CO v(n,P, r s , k) is a piecewise constant function ofr s . In particular, there 
exist = ro < r\ < r-i < ... < r m < -j=, and c\, C2, c rn such that p cov (n,p 7 r s ,k) = c-i if 

r s e [ri-!,ri). 

PROOF. For simplicity we prove the theorem for k = 1; the extension to k > 1 is 
straightforward. Let p CO v(r s ) = p CO v(n,p, r s , 1). We need to show p co v(r s ) ia a piecewise 
constant function of r s . It is clear that p CO v(r s ) is a nondecreasing function. In partic- 
ular we have p CO v(0) = and p CO v(r s ) — 1 for r s > -±=. For a point X in the plane, let 
circ(X, r) be the closed ball that is centered at X and has radius r. Define cov(X, r s ) 
to be the area that is covered by a sensor node located at X with sensing radius r s . 
In other words, cov(X, r s ) is the portion of circ{X, r s ) that lies within the deployment 
region. Again assume that X = {Xi,X 2 , ...,-Xjfe} is the set of points in the sensor grid 
and X a c X is the set of active sensors. Define 

cov(X a ,r s )= (J cov(X,r s ). (3) 
xex a 

Thus the unit square So is completely covered whenever cov(X a ,r s ) = So- Let < r s < 
If cov(X a , r s ) = Sq, then for all r > r s , we have cov(X a ,r) = S . On the other hand, 

we prove that if cov(X a ,r s ) ^ S , there exists e > such that for all r e [r s ,r s + e) we 
have cov{X a , r) ^ S Q . To prove this note that the covered area cov{X a ,r s ) is a closed set 
because it is the union of a finite number of closed sets. Thus, the uncovered area is an 
open set. Hence, to cover the uncovered area, the sensing radius r s must increase by a 
strictly positive amount. 

We now prove that for any r s , there exists a strictly positive e such that p cov (r) 
remains constant as the sensing radius r varies within [r s ,r s + e). Define 

X r a ° = {X a : cov{X a ,r s ) = 5 }, and = {X a : cov(X a , r s ) ^ S }. (4) 

Note that X^ s and X^ s are finite sets. Using (HJl we have 

Pcov{r s ) = ^2 p(X a ). (5) 

For any X a e X a s define e{X a ) = min{r' s.th. cov(X a ,r s + r') = S } and let e = 
min{e(A" a ) : X a <s X^ s }. Then e > 0. Further, for all r e [r s ,r s + e), we have X^ = X^\ 
Thus we conclude that for all r e [r s ,r s + e), we have X£ = X^ s . Using (O we conclude 
that p CO v{r s ) does not change as r varies in [r s ,r s + e). This proves that p CO v(r s ) is a 
right-continuous piecewise constant function. 

It remains to show that the number of discontinuities is finite. This follows easily 
from the fact that the number of X a 's is finite. Note that by (©, any discontinuity occurs 
when the set X£ s changes due to an increase in r s . However, X£ s can have at most 2" 
elements. Further, at each discontinuity, at least one element is added to XJ/. This 
implies that the number of discontinuities is upper-bounded by 2". It is worth noting 
that in practice, the number of discontinuities is much smaller than 2™. This completes 
the proof. □ 

Theorems|2]and[3]determine the behavior of a vast class of network quantities when 
they are considered as functions of communication and sensing radii. In particular, 
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these are important from the view point of finite sensor grids. We note that for very 
large network sizes, the piecewise constant functions tend to continuous functions. 
Thus, we do not observe the discontinuities. However, in such networks as finite sen- 
sor grids, this property is noticeable as in Figures M and [3] We clarify that because 
the simulation results are approximations for the actual values, the figures are not 
completely piecewise constant. In fact one of the implications of Theorems [2] and [3] is 
to simplify simulations since the piecewise constant functions can be completely de- 
termined by knowing their values for only a finite number of points. Furthermore, 
the above results suggest that increasing the communication and sensing radii does 
not necessarily improve coverage, connectivity or any other graph theoretic properties. 
This is an important observation for designing the network and choosing its parame- 
ters optimally. 



2.3. Bounds on the Coverage Probability 

We now consider coverage probability for finite sensor grids. We find lower and upper 
bounds for p cov (n,p,r s , k) and show that they can give an acceptable estimate of the 
coverage probability. Let N(r, x, y) be the number of sensors whose distance from the 
point (a;, y) is less than or equal to r. For example, N(r, 0.5, 0.5) denotes the number of 
sensors whose distance from the top-right corner of unit square is less than or equal 
to r. We first prove the following lemma. 



LEMMA 1. Let L be the set of ' \fl x \fl points in a virtual grid on the unit square. 
Let us also denote by A{u) the event that point u on L is covered by a sensor grid with 
coverage radius r s . We then have 



Proof. We use FKG inequality to prove this lemma MFortuin et al. 197 ll . We first 
show that for any two subsets I and J of L, we have Pr(A(I) f\ A(J)) > Pr(A(I)) x 
Pr(A(J)). Here, A{I) (A(J)) is the event that all points in / (J) are covered. Since this 
is true for any two subsets of L, © can be derived by partitioning L and the resulted 
components, repeatedly, and then using this property at each step. 

First note that we can enumerate the nodes in the sensor grid from 1 to n. Accord- 
ingly we can show the status of the network with a n-tuple binary vector where and 
1 are assigned to inactive and active nodes, respectively. Let us denote by T the set of 
all possible binary ?i-tuples as the network status, i.e. T = {t = (ti, t 2 , t n ) e {0, 1}"}. 
T can be then defined as a "finite distributive lattice" as follows. For x = (x\, x 2 , x n ) 
and y = (j/i, y 2 , —, y n ) in T, we define x V y as the elementwise "or" of x and y. That 
is if w — x V y = (wx,w 2 , w n ) then w 8 = i, V yi. Similarly, we define x A y as the 
elementwise "and" of x and y, i.e. if w = x A y = (wi,w 2 , w n ) then w, = i, A yi. With 
these definitions, it is easy to check that V and A are distributive over each other. Note 
that the lattice defined this way is partially ordered as we have (x A y) ■< x, y -< (x V y). 

We now define a probability measure fx : T -> R + as follows. For x e T, fi(x) = 
p k (l - p) n ~ k , where k = Yn=i Xi - Note that p(x.) in fact indicates the probability that 
the sensor grid admits the status x with k active sensors and n — k inactive sensors. 
It is also trivial to verify that fi(x)p,(y) < /i(x V y)^(x A y), which is required by FKG 
inequality. Given two subsets I and J of L, we also define functions /, g : T — > R+ as 
follows. For every x e T, /(x) = 1 (g(x) = 1) if I (J) is covered. By these definitions, 
/ and g are both increasing functions over T. Given the lattice T, measure \i, and 




(6) 
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functions / and g as above, the FKG inequality holds as follows. 

(]>>(x)/(x)).(5>(x).9(x)) < (5>(x)/( X ) 5 ( X )).(5>( X )). (7) 
xgt xgt xgt xgt 

However, £ x6T M x ) = !• Furthermore, E xS t M x )/( x )> ExeT A*( x )5( x ), and 
£ xeT /*(x)/(x)g(x) are in fact equal to Pr(A(/)), Pr(A(J)), and Pr(A(/)y/\A(J)), re- 
spectively. © can thus be rewritten as Pr(A(I) f\ A(J)) > Pr(A(I)) x Pr(A(J)). This 
completes the proof. □ 

Now, we are ready to prove the lower bound on the coverage probability. 

THEOREM 4. Consider the coverage probability for a finite sensor grid with param- 
eters n,p, and r s . We then have 

Pcov (n,p,r a ,l) > l[[l-(l-p) N ^>^], (8) 

where L is the virtual grid in Lemma\l\ and the radius r' s is given by r s — -^=. 

PROOF. First note that the choice of the virtual grid L and its size, I, is arbitrary. 
As a result, for any given r s , we choose / large en ough such that r s — \j\f% > 0. To 
prove this theorem, we make use of some results in Kumar et al. 20081 . Lemma 3.1 in 
[Kumar et al. 20081 states that for a given set of points L that consists of all grid points 
of a Vl x \fl virtual grid on a unite square, if L is covered by a network of radius r' s , the 
unit square is covered by the same network but with the radius r s = r' s + Hence, 

Pcov{n.,p, r s , 1) > Pr(L covered). Now we use Lemma[T]above to prove the lower bound. 
Let us denote by A{u) the event that point u is covered by a network with coverage 
radius r' s . By Lemma Q] we have Pr(/\ uei A(u)) > YIuel Pr(A(u)). Now, note that the 
probability that a point u with coordination (x u , y u ) is covered by the set of n nodes with 
coverage radius r' s is given by [1 - (1 - p) N i r ^ x ^'V^)]. Thus, we write p CO v(n,p, r s , 1) > 
Pr(L covered) = Pr(A„ eL A(u)) > U ueL Pr(A(«)) = ft^l - (1 - p)^:.*«.w-)]. n 

Now, we prove an upper bound for the coverage probability. 

THEOREM 5. Consider sensor grids with parameters n, p, r s . Then the coverage 
probability is upper bounded by 

p cov (n, P , r s , 1) < [1 - (1 - p)^('"-8,5) ] 4 X [1 _ (1 _ p) iV(,« ) . S ,0) ] 4 L ll^ij 

x [l-(l-p) Ar ^^°)]L ii 2^ 1 J 2 , (9) 

where [x\ denotes the largest integer less than or equal to x. 

PROOF. LetXi = (xi,j/i),X 2 = (x2,V2),—,X m = (x m ,y m ) be m points on the deploy- 
ment region So. Assume that d(Xi,Xj) > 2r s for i ^ j, where d(., .) is the Euclidean 
distance between the points. Then the event that Xi is covered is independent of the 
event that Xj is covered. This is because there is no sensor node that can cover both 
points. Hence, the probability that all X/s are covered is given by 

m 

- (1 - p)^-.»«*)]. (10) 

i=l 

This implies that p cov (n,p,r s , 1) is upper bounded by n2=i[l ~ (1 — p) Ar ( r s,a:i,j/»)] - Thus 
using any set of points on the plane that satisfy Xj) > 2r s , we can find an upper 
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bound for p cov (n,p, r s , 1). In particular, considering the set of points given by Figured! 
we obtain the upper bound in ©. □ 




Fig. 4. Location of the points used for the upper bound of theorem [5] The centers of the circles are points 

Xi = {xi,yi),X 2 = (x 2 ,V2), —,X k = {x k ,y k ). 

Note that the choice of Xi's in the proof ensures that we consider the edge effects. In 
fact, in many situations, the coverage probability is dominated by the first and second 
terms in © which are related to edge effects. One may suggest that using a triangular 
grid, instead of non-overlapping balls, can result in a more dense packing and conse- 
quently a better bound. However, using a triangular grid results in fewer nodes on the 
sides of the square. We have evaluated dlOD for the triangular grid as well as some 
other more complicated layouts. It turned out that the resulting bound is looser for 
the triangular grid. Moreover, there is only a negligible improvement by using other 
layouts at the expense of a more complicated expression compared to ®. It is also 
worth noting that N(r s , .5, .5), iV(r s , .5, 0), and 7V(r s ,0, 0) introduce discontinuities in 
the upper bound as predicted by Theorem [3] 

Figure [5] compares the results obtained by Theorems |4] and [5] and the simulations 
for n = 100 and p = 0.2. We observe that Theorems |4] and [5] provide significantly better 
estimates of coverage probability compared with the asymptotic analysis in Figure [2j 
The asymptotic behavior of these bounds can be checked by letting n to grow large. The 
de rivation of the lower bound employs a similar argument as in the case of Lemma 4.1 
in [Kumar et al. 20081 . It can be checked that this bound leads to the same asymptotic 
expression as in Theorem [H hence it is tight asymptotically. On the other hand, when 
n gets large, we can reasonably expect the same situation as in the upper bound of 
Theorem[5j That is the terms corresponding to the virtual nodes (nodes on the virtual 
grid) on the corner and close to the edges will be dominant in the lower bound of (©. 
This is true because there are fewer sensor nodes around these virtual nodes to cover 
them, causing the coverage probability Pr(A(u)) for these virtual nodes to decay faster 
than the rest of the virtual nodes. Regarding the asymptotic behavior of the upper 

2 

bound of (O, it can be verified that if iog(n P ) < 1 ~ e as n tends to infinity, then the 
upper bound will be o(l). 
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Sensing Radius r s (n=100, p=0.2, k=1) 

Fig. 5. Comparison of finite-size analysis with the simulation results for the coverage probability of sensor 
grids with n = 100 and p = 0.20. 



We also like to talk about the time complexity of computing the bounds in Theorems 
|4]and[H The u PP er bound of Theorem[5]can be computed in time 0(n). This is because 
we need to find the neighbors for l/(2r s ) 2 points, and finding the number of neighbors 
for each point takes a constant amount of time. However, if 1/ (2r s ) 2 > n, then p cov = 0. 
Thus, the complexity is 0(l/(2r s ) 2 ) = 0{n). For the lower bound in Theorem |H note 
that we only need to find the number of neighbors for every node of L. Given the sensor 
grid and the virtual grid L, finding the number of neighbors for each node of L takes a 
constant amount of time. Also note that for r' s to be positive, L needs to contain more 
nodes than the sensor grid, hence, I > n. Therefore, the lower bound can be computed 
with complexity 0(1). 

Theorems |4] and [5] can be easily generalized for fc-coverage. Since the proof is very 
similar, we just state the result in one theorem. 



THEOREM 6. Consider the k-coverage probability for a sensor grid with parameters 
n, p, and r s , and assume that L and r' s are as defined in Theorem^ Then we have 

p cov {n lP ,r s ,k)>X{[l- k f^( N{r ^^ (11) 
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- Simulation Results 

■ Lower Bound 




0.4 0.6 0.8 1 

Sensing Radius r (n=100, p=0.2, k=2) 



Fig. 6. Comparison of finite-size analysis with the simulation results for the 2-coverage probability (k = 2) 
of sensor grids with n = 100 and p = 0.20. 



and 



fc-1 



,(n,p,r s ,k) <[1 - ^2 



i=0 
fe-1 



fN{r s , .5, .5) 



^(l-p) Ar ( , '- 5 <- 5 )- l ] 4 x 



[1-J2 ( N( - rs v 5,0) ) P i (i - P ) N{r - 5 ' 0) ^' llL1 ^ li 



i=0 
fc-1 



i=0 



/iV(r s ,0,0) 



(12) 



Figure [6] compares the results obtained by Theorem[6]and the simulations for k = 2, 
n = 100, and p = 0.2. We observe that the two results are very close. 

3. SMALL-SCALE ANALYSIS FOR RANDOM SENSOR NETWORKS 

In this section, we try to establish a framework for analysis of finite sensor networks 
with random node deployment. As we mentioned earlier, the exact analysis of network 
properties is usually very difficult or at least results in very complicated formulas. 
Thus, we will try to find simple lower and upper bounds which are sufficiently close 
together that can be used to find a good estimate of the exact value of the desired 
property. Here, we consider coverage and connectivity in finite sensor networks. 

3.1. Preliminaries 

We consider a wireless sensor network that consists of n nodes and assume that the 
nodes are placed on a plane based on a given probability distribution. For example, in 
wireless sensor networks it is usu ally assumed that th e nodes are randomly and uni- 
formly deployed over a given field | Akyildiz et al. 2002 1. We assume that each node has 
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1 



r(n) 







S3 















Fig. 7. The field So and its subdivisions 81,82, and 53. 

a fixed communication radius. Two nodes are connected (can communicate with each 
other) if they are within communication range of each other. Throughout the paper, we 
assume B(M?) is the Borel er— algebra on R 2 and m is the Lebesgue measure on S(M 2 ). 
Note that we just use measure theoretic definitions to take care of technicalities but it 
is not necessary for the reader to be familiar with them. The reader can simply assume 
that for a set F in R 2 , m(F) is the area of F. B(X, R) is the closed ball with radius R 
centered at X in R 2 . S(X, L) is the closed square with side L centered at X in R 2 . In 
particular So = S(0, 1) is the closed square with unit area centered at the origin. If u 
and v are two nodes of a network located in M 2 , then d(u, v) is the Euclidean distance 
between the location of the points. For any set F e B(R 2 ) we define v{F) = m(F n So). 
Clearly, v defines a measure on $(R 2 ). 

Wireless networks are sometimes modeled with the probability space of graphs that 
we represent with g(n,r) = g{n,r{n)). In this model, it is assumed that n nodes are 
uniformly and randomly distributed over So = S(0, 1). If two nodes u and v satisfy 
d(u, v) < r(n), then the edge {u, v} belongs to edges of the graph. A more general model 
is the model g(n,r,p), in which two nodes are connected with probability < p ^ 1 
if their distance is less than r. In this model p models link failures that are common 
in wireless networks. Note that here we are using p as a different notation from the 
previous section. Asymptotic properties of g(n,r) have been studied extensively Here 
we are interested in these properties when n is not necessarily large. It is worth noting 
that the assumption that the nodes are distributed on a square is made for simplicity 
These arguments can easily be generalized to other models for the deployment region 
as well as the case where nodes are distributed non-uniformly over the deployment 
region. For the purpose of analysis, we divide the square So to different parts shown in 
Figure El 

3.2. Asymptotic versus Finite Analysis 

In this section, we present some evidence to show that previous asymptotic results 
diverge significantly from actual values for finite networks. To show this, we con- 
sider connectivity We first provide the asymptotic probability of disconnectivity for 
g(n,r,p) and compare it to simulation results. The following result is proved in 
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Fig. 8. Comparison of asymptotic results with the small scale simulation results for the probability of 
disconnectivity of g(n = 100, r). 



[Gupta and Kumar 1998], where a slightly different model is considered. However, the 



results can be trivially extended to g(n, r). 

THEOREM 7. (Gupta and Kumar 1998) Let c n = mxr 2 - log(n), then g(n,r) is con- 
nected with high probability if lim c„ = oo. On the other hand, if lim c„ = c < oo then 

for large n, g(n, r) is disconnected with a strictly positive probability 1 — p aS ymp{c). 

This theorem states that if lim c„ = c < oo, the network connectivity probabil- 

ity will be bounded away from one. In fact, p a symp{c) is the limit for the probability 
that the netwo rk is connected when n goes to infinity. To find p a sym P (c), Penrose in 
[Penrose 19971 proved that g(n,r) is connected if and only if the longest edge of its 
corresponding Minimal Spanning Tree (MST) is smaller than r. On the othe r hand, if 
we denote the longest edge of the MST by M n , it is shown in BPenrose 1997B that the 
distribution of mrM n 2 — logn converges to the double exponential distribution: 

lim P[rmM n 2 - logn < a] = exp{-e~ a ) for a G K. (13) 

n— ¥cg 

Thus we have 

Pas ymp (c) = lim P[M n <r]= lim P[mrM n 2 - logn < nirr 2 - logn] = e _e_c . (14) 

n— >oo n— foo 

Therefore, asymptotically, the probability that g(n, r) is connected is given by 

n — e ~ ne ~"*"' 2 

Fasymp c 

In Figure [H we compare the probability of having a disconnected graph for n = 100 
and for both exhaustive simulations and the asymptotic results. In Figured! the prob- 
ability of disconnectivity is shown as a function of r, the communication radius. The 
experiment shows that these results may differ by 10 orders of magnitude. This il- 
lustrates that the asymptotic method fails to provide a good approximation for small- 
scale networks. 
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A natural question to ask is what makes the results for the asymptotic analysis so 
different from the finite case? As you can see in Figure (7) So is formed by three re- 
gions, Si, and boundary regions 5 2 and 6*3. One important phenomenon in asymptotic 
analysis is that boundary effects can be neglected. Loosely speaking, the asymptotic 
analysis of the network properties is usually dominated by what happens in region 
Si in Figure [7] This can considerably simplify the analysis and results in simple and 
closed-form formulas for network properties. In fact, we saw an example of this phe- 
nomenon in the asymptotic formula for connectivity in d!4l >. However, in small-scale 
networks boundary effects cannot be neglected. In other words, nodes in the corners of 
the field can play an important role in some network properties. 

Another important issue in the analysis of finite networks is the effect of constant 
factors. In asymptotic analysis, we usually neglect constant factors. However, in the 
small-scale analysis, we must consider them. This is in fact a distinction of any finite 
analysis from the asymptotic analysis and is not specific to geometric graphs. 

3.3. Small-Scale Analysis for Coverage 

In this section, we study the coverage probability, p cov (n, r), for finite sensor networks 
modeled by g(n,r). We prove lower and upper bounds for the coverage probability. 
We start with the lower bound which gives the worst case performance as well as a 
guarantee of the coverage probability. 

THEOREM 8. Consider the coverage probability of a sensor network modeled by 
g(n, r). Then we have 

Pcov(n, r) > 1 - £[1 - v(B(X(u),r'))] n , (15) 

where L is the set ofy/lxy/l points in a virtual grid on the unit square and the radius 
r' is given by r — -^=. 

Proof. We briefly describe the proof. As in the case of grid deployment, showing 
that L is covered guarantees the coverage of the entire region. Let A(u) be the event 
that the virtual grid point u is covered. Using union bound, we have 

Pr(L not covered) = Pr[ (J A(uj] < ^ Pr(A(^)) = ^[1 - u(B(X(u),r'))] n . (16) 

Therefore, p cov {n,r) = 1 - Pr(L not covered) > 1 - £ ueL [l - v( B ( x ( u ), r 'W- □ 
Now we prove an upper bound for the coverage probability. 

THEOREM 9. The coverage probability of a unit square for a sensor network mod- 
eled by g(n, r) has an upper bound given by 

2 2 

Pcov(n,r) < [1 - (1 - ^H 4 x [1 - (1 - ^-)™] 4 L^J x [1 - (1 - irr 2 )"}^^^ . (17) 

Proof. We adapt the proof of Theorem [5] to prove Theorem [9l Consider k points 
Ui, U2, Uk on the unit square and assume that these k points are at least apart by 
2r units from one another. Similar to the proof of Theorem [5) we can observe that 
Pcov(n, r) is upper bounded by the probability that all the k points are covered which is 
given by n!Li(l ~ [1 ~ v {B{Ui, r))] n ). Using the set of points depicted in Figure |H we 
find the upper bound given by dlYI . □ 

Figure [9] compares the bounds predicted by Theorems [ 8] and [9] with th e simulated 
coverage probability value. Asymptotic result from [Kumar et al. 20081 is also pre- 
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Fig. 9. Simulation results, and upper and lower bounds of coverage probability for a random wireless net- 
work of size n = 100. 



sented. Clearly, the bounds are more useful than the asymptotic result in the sense 
that they give a better estimate of the coverage probability. 

3.4. Small-Scale Analysis for Connectivity 

In this section, we study the connectivity properties of finite sensor networks mod- 
eled by g(n, r,p). We find lower and upper bounds for the probability Pdisc(n, r,p) that 
g(n, r, p) is disconnected. Let pi ow (n, r, p) and p upp (n, r, p) be the lower and upper bounds 
on Pdisdn, r, p), respectively. Here we consider the case where Pdisc(n, r, p) is small, i.e., 
Pdisc(n,r,p) < 0.1. In practice, this is usually the range that is important, since we 
want the network to be connected with high enough probability. Using these bounds, 
we then provide a simple formula to estimate pdt SC (n, r,p). As we will see by simula- 
tions, the proposed formula gives a very good estimate for pdi SC {n, r,p). First, note that 
a connected component of a graph g is defined as a connected subgraph that is isolated 
from the rest of g. 

THEOREM 10. Consider a wireless sensor network modeled by g(n,r,p). Then we 
have 

Pdisc(n,r,p) >n J ^1 - v{B{X,r))p S j dm(X) 

So 

So So 

?i-2 

+ u(B(X,r) nB(Y,r))p 2 ) dm(X) x m(Y), (18) 
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and 

n/2 



Pdisc(n,r,p) (j^jpcomp({vi,v 2 , —,v k }) =n J ^1 - v{B(X, r))p^ dm(X) 



k=1 So 
n/2 

' n 



k=2 



^2 ( h }Pcomp{{V\,V2, -,V k }), (19) 



where p C om P ({vi,v 2 , v k }) is the probability that the vertices in {vi, v 2 , v k } construct 
a connected component in g(n, r,p). 

Proof. Let pi(n, r,p) be the probability that there exists at least one isolated node 
(a vertex with no neighbors) in g(n,r,p). Let also Vi,V2,.-,v„ be the n vertices of 
g(n,r,p). Then Pdisc(n,r,p) > pi(n,r,p). Applying the inclusion-exclusion lemma we 
obtain 

Pdisc(n,r,p) >y^(-l) fc+1 ( J VY{vi,v2,...,v k are isolated vertices} 

fc=i ^ ' 

>n Pr{vi is isolated} - ^^Pr{ui and v 2 are isolated}. 

Note that Pr{vi is isolated } = / ^1 - u(B(X,r))pj dm(X). Now define 
Circ(a, b, r) — {(x, y) : (x — a) 2 + (y - b) 2 < r 2 }. Then we have 

Pr{«i and v 2 are isolated vertices} = J J ^1 — v(B{X, r))p - v(B(Y, r))p + 

Sa S \Circ(X,r) 

p{B(X, r) n B(Y, r))p 2 ) dm(X) x m(Y) + (1 - p) x 



i n-2 
2 ' 



1 - v(B{X, r))p - u(B(Y, r))p + v{B(X, r) n B(Y, r))p z ) dm{X) x m(Y) 

Sa Circ{X,r) 

t n-2 
2 ^ 



So So 

Combining these equations, we conclude the lower bound. For the upper bound, note 
that pdisc{n, r,p) is equal to the probability that g(n, r,p) has at least one component of 
size less than n/2. This is given by eq. ( fl9l . □ 

Note that the bounds for Pdisc{n,r,p) may not satisfy the simplicity requirement. 
Particularly in the upper bound, except the first few terms, finding the rest of them is 
computationally infeasible. We now try to give an estimation of pdisc(n,r,p) based on 
these bounds. Let us denote the fcth term in the upper bound by a k . We recall the as- 
sumption that pdiscin, r,p) is not very large, specifically we assumed Pdi sc (n, r,p) < 0.1. 
An important observation here is that, by this assumption, the a k coefficients decay 
very fast and, hence, the term J2 k =i a k is dominated by a\. This can be seen by both 
numerical sim ulations and intuitive analytical arguments. In fact, as it is shown in 
Penrose 2003L as n tends to infinity, the impact of the terms a k , k > 1 fades. Figure [101 
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Fig. 10. Comparison of a± and ao in H9\ . 

compares a\ and a 2 for g(n — 100, r,p = 0.5). As we see a 2 is at least one order of mag- 
nitude smaller than oi. Using the same approach, we find out that a similar argument 
is true about the first and second terms in the lower bound of eq. ( fl8l l. However, the 
first term is shared by both the lower and upper bounds. Based on these observations 
we approximate the probability of disconnectivity as follows. 



Pdtsc(n, r,p) ~ n / 1 - v{B(X, r))p dm(X) 



(20) 



Figure [TT] shows the upper bound, lower bound, and the simulation result for the prob- 
ability of disconnectivity of g(n, r,p), for n = 100, and p = 0.5. As it can be seen, the 
three curves almost overlap. Based on our simulations, similar results are achieved if 
we use different choices of parameters. 

It is worth noting that the methodology used here can be used to study k- 
connectivity. In summary, we find the following approximation of the probability that 
g(n, r) is not fc-connected 



p k , disc (n, r) ~ n (") / M W x ( 1 - V W X ' 

•So 



1=0 



n—j—1 



dm(X). (21) 



Our simulations for different values of k confirm the validity of d21l >. Here, due to the 
space limitations we omit those results. 

4. CONCLUSIONS 

In this paper, we took some initial steps towards analyzing finite wireless sensor net- 
works. We provided some compelling evidence to show that asymptotic results are not 
suitable for analyzing practical finite sensor networks. We studied connectivity and 
coverage of finite unreliable sensor grids as a special case. We showed that the connec- 
tivity, as well as all the graph theoretic quantities, are piecewise constant functions 
of the transmission radius in such networks. We also proved that the coverage has a 



ACM Transactions on Sensor Networks, Vol. 10, No. 1, Article 0, Publication date: 2014. 



0:20 



A. Eslami et al. 




Fig. 1 1 . Disconnectivity probability of g(100, r, .5): lower bound, upper bound, and the simulation results. 



similar behavior. Moreover, we obtained lower and upper bounds for the coverage and 
fc-coverage probability of the grids and verified their preciseness through simulations. 
Next, we extended our study to finite sensor networks with random node deployment. 
Specifically, we considered coverage and connectivity of such networks. We derived 
lower and upper bounds for their coverage and showed how they can be used to esti- 
mate the coverage probability of the network. We also obtained a formula for connec- 
tivity of wireless sensor networks and verified its accuracy through simulations. The 
formula was then extended to include fc-connectivity. A common characteristic of all 
these bounds is the ease of computations, making them very attractive. 

This paper also opens up many research possibilities that offer some potentials for 
further study. In the past, many other important properties of wireless sensor net- 
works have been studied for large-scale networks. It is an important task to extend 
these results for networks with practical sizes, i.e. small-scale networks. Small-scale 
analysis can also reveal the effects of network parameters on network characteristics. 
The next step would be to derive more accurate bounds for network parameters such as 
coverage, connectivity, and MAC layer capacity and further use the small-scale frame- 
work in the design, analysis, and evaluation of communication algorithms for wireless 
networks. 
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